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INFLUENCE  OF  SHEAR  MODULUS  ON  THB  BEHAVIOR  OF 
ACOUSTICALLY  TRANSPARENT  MATERIALS 

INTRODUCTION 

Several  publications  [1-4]  have  appeared  in  recent  years  that  report  on 
the  development  of  an  acoustically  transparent  material  based  on  a  fluoroepoxy 
filled  with  microballoons.  Such  a  material  can  be  designed  to  have  a  density 
and  dilatational  sound  speed  that  closely  match  those  of  the  acoustic  medium, 
usually  seawater,  but  with  sufficient  rigidity  to  serve  as  a  material  for 
structural  members  or  backing  plates  of  acoustical  systems.  The  matching  of 
density  and  dilitational  sound  speed  leads  to  perfect  acoustic  transparency 
for  reflection  of  sound  at  the  interface  of  fluids  only.  A  solid  has  a  finite 
shear  modulus;  and  for  other  than  perpendicular  incidence,  a  shear  wave  will 
be  present  in  addition  to  the  dilatational  wave.  This  shear  wave  causes  a 
deviation  from  ideal  transparency. 

In  this  study  an  analysis  is  given  of  the  influence  of  the  shear  modulus 
on  the  reflection  and  transmission  of  a  plane  sound  wave  by  an  infinite  plate 
with  strictly  or  approximately  matching  density  and  sound  speed.  It  shows  the 
limits  on  plate  thickness  and  frequency  imposed  by  the  requirement  of 
negligible  reflection.  The  effect  of  a  finite  loss  factor  in  the  shear 
modulus  is  also  investigated. 

ANALYSIS 

The  analysis  of  the  reflection  and  transmission  of  a  plane  sound  wave  by  a 
plate  is  given  in  Ref.  5.  The  main  points  of  the  derivation  are  repeated 
here.  In  the  general  case  (Fig.  1)  a  fluid  with  density  pQ  is  present  at  one 
side  of  the  plate  (thickness  h  =  2d)  where  the  incident  wave  arrives,  with 
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Fig.  1  -  CcoHtrj  of  f luid-lo»d«d  plots 


angle  8  to  the  normal  and  wavenumber  kQ.  Another  fluid  with  density  p'  and 
sound  speed  c1  is  at  the  other  boundary  of  the  plate.  The  transmitted  wave 
forms  an  angle  8 '  with  the  normal. 

The  propagation  of  waves  in  the  solid,  excited  by  the  incoming  wave,  is 
treated  by  exact  elasticity  theory.  The  plane-crested  wave  in  the  solid  with 
wavenumber  k  is  a  combination  of  a  dilatational  wave,  represented  by  a  scalar 
potential  ^  with  amplitudes  A  and  B  for  the  symmetric  and  antisymmetric 
contributions,  respectively,  and  a  shear  wave  represented  by  the  s-component 
of  a  vector  potential  f  with  amplitudes  D  and  C  .  Other  field  variables  are 
the  pressure  amplitudes  of  the  incoming  wave  P^,  of  the  transmitted  wave  P' , 
and  of  the  reflected  wave  Pf.  At  the  boundaries  of  the  plate  with  the  two 
fluids,  one  imposes  conditions  of  continuous  normal  stress  and  velocity  and 
sero  tangential  stress.  The  matrix  of  the  coefficients  of  the  seven  field 
variables  in  the  six  boundary  conditions  (after  manipulation  of  the  rows)  is 
shown  in  Table  1.  Here  P  =  P.  exp  (ik  d  cos  6) ,  P  =  P  exp  (-ik  d  cos  8), 


Matrix  of  Coefficients  of  Equations  Describing  Reflection  and  Transmission  by  Fluid-Loaded  Plate 


P,*,IT"Pjn  PibaI*P,b 
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and  P'  =  P.  exp  (-ik'd  cos  8]).  The  wavenumbers  k  ,  k,  and  k'  are  related 

—  ®  2  2  2 

through  the  coincidence  condition  k  =  k  sin  6  =  k'  sin  6* .  Also  q'  =  k ,-k 

2  2  2  ®  ® 
and  s'  =  k  -k  .  where  k,  and  k  are  the  wavenumbers  of  the  dilatational  wave 
s’  d  s 

and  shear  wave  in  the  solid. 

The  ratios  of  the  field  variables  may  be  found  from  the  matrix  as  ratios 
of  the  corresponding  subdeterminants.  One  finds  for  the  case  where  the  same 
fluid  is  present  on  both  faces  of  the  plate  that  the  reflection  and 
transmission  coefficients  R  and  T  are  given  by 


(2ikQd  cos  6) 


a: 


T 


exp  (2ik  d  cos  8) , 

«  ■?  *  r\ 


(2: 


where 


2 _ -  A  +  D 

p  ~A  +  B-  C  +  D’ 
o 


(3: 


sL  =  b  +  c 

p  A  +  B  -C  ♦  D  ’ 
o 


(4: 


*4 


2 


2  ^ 

(q' d)  (k  d)  cos  q' d  sin  q1  d  cos  s'd  sin  s'd, 


(5: 


k  d  q'd  t 


cos  q’ d  cos  s'd  cos  8, 


(6' 
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C  = 


k  d  q'd  A 
o  n  a 


sin  q'd  sin  s'd  cos  6, 


(7) 


D  =  i  (kod) 2  4a  4s  (cos2  «)/(ksd)4, 


(8) 


and 


Aa  =  “ll"^  "  m12m21  ’ 


(9) 


As  =  m33m44  "  m34m43  ’ 


(10) 


in  terms  of  the  elements  m. .  of  the  matrix  in  Table  1. 

ij 

It  is  of  interest  to  formulate  the  condition  for  zero  transmission  through 
the  plate,  which  is  equivalent  to  total  reflection  if  there  are  no  losses  in 
the  plate.  From  Eq.  (4)  one  sees  that  this  is  equivalent  to  B  +  C  =  0. 
Inserting  the  expressions  for  A  and  A  from  Eqs.  (9)  and  (10),  one  may  reduce 
this  to 


[(kd)2  -  (s'd)2]2 
♦  4(kd)2 


cos  s'd  sin  s'd 

q'd  s'd  cos  q'd  sin  q'd  =  0. 


(11) 


A  computer  program  was  written  to  evaluate  the  reflection  coefficient  R 
and  transmission  coefficient  T  according  to  Eqs.  (1)  and  (2).  Input  to  the 
program  is  the  ratio  of  density  of  the  fluid  to  that  of  solid  P0/Ps>  sound 
speed  in  the  fluid  relative  to  the  dilatational  speed  in  the  solid  co/c^,  the 
ratio  of  shear  to  bulk  modulus  in  the  solid  G/K,  the  loss  tangent  in  the  shear 
modulus  f?Q,  and  the  thickness  of  the  plate  relative  to  the  wavelength  of 
dilatational  waves  in  the  plate  material  h/X.  The  bulk  modulus  K  is  assumed 
lossless . 

If  one  introduces  the  following  symbols 
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2 

a  =(/>//>)  q'd  (k  d)  cos  q' d  cos  s'd 
r  o  r  s  s 

a'  =  (/>0//>g)  q' d  (kgd)2  sin  q' d  sin  s'd 

d  =  A  k  d/(k  d)2 

d'  =  A  k  d/(k  d)2  , 
s  o  '  v  s  / 


one  may  write  for  the  expressions  A,  B,  C,  D, 

A  =  1/2  aa' 

B  =  1/2  i  ad'  cos  0 
C  =  1/2  i  a'd  cos  8 
D  =  1/2  dd1  cos2  8  . 


The  reflection  factor  P/Pq  may  be  expressed  as 


P_  =  (d'/a'Ud/al  cos*  8  -  1 
Po  (1  -  i  —  cos  8) (1  +  i  ^7  cos 


A  "structural  response  function"  ft  was  introduced  in  Ref.  6 


ill  cos  8  - 


P  +  P 
o 

P  -  P  • 
o 


by  setting 


Here  one  finds  a  similar  structural  response  function  from 
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o  _  (d1  /a')  -  (d/al  -  2ifd/aHd'/al)  cos 
2  +  i (d*  /a.' )  cos  8  -  i(d/a)  cos  8 


The  structural  response  function  is  related  to  the  effective  impedance  Z 
of  the  plate  surface,  defined  as  the  ratio  of  the  total  pressure  P  +  Pq  at  the 
side  of  the  incoming  wave  to  the  surface  particle  speed  perpendicular  to  the 
plate.  One  finds 


pc  P  +  P 

rj  rO  O  O  .  n  „ 

2  "  '  cos  8  P  -  P  "  -l(1W 
o 


One  may  define  effective  surface  impedances  for  antisymmetric  and 
symmetric  waves  separately.  First  the  pressure  at  the  side  of  the  incoming 
wave  is  split  in  a  part  contributing  to  antisymmetric  waves  and  a  part 
contributing  to  symmetric  waves,  namely 


=  (p*F0*r'). 


Similarly  the  displacement  perpendicular  to  the  plate,  w(d) ,  is  written  as 


*(d)  =  \  [w(d)  ♦  w(-d)]  ♦  \  [w(d)  -  w(-d)]  . 


The  two  impedances  are  then  defined  by 


Pq  ♦  F  -  P' 
2a  _  w(d)  +  w(-d) 
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Z  = 


P  +  P  +  P' 

o  _ 


s  w(d)  -  w(-d) 


(21) 


In  terms  of  the  symbols  introduced  before,  one  finds  Z&  =  -id/a  and 

Z  =  id' /a'  . 

s  ' 


It  would  be  instructive  to  express  Z  in  terms  of  Z  and  Z  in  the  form  of 

r  as 

an  equivalent  circuit  diagram.  Thus  far  this  has  not  been  successfully 
accomplished. 


RESULTS  AND  DISCUSSION 

Figure  2  shows  the  reflection  coefficient  as  a  function  of  incidence  angle 
6  for  various  values  of  relative  density  and  sound  speed.  The  three  figures 
arranged  along  a  column  have  equal  values  of  p  /p  ,  namely  0.995,  1.000,  and 
1.005.  The  three  figures  arranged  in  a  row  have  equal  value  of  co/c^:  0.995, 
1.000,  and  1.005.  The  four  curves  marked  1,  2,  3,  4  vary  in  value  of  the 
thickness  h  (=2d)  relative  to  the  wave  length  X  of  dilatational  waves  in  the 
material.  The  four  curves  correspond  to  h/X  =  0.50,  0.25,  0.10,  0.05.  The 
value  of  the  shear  modulus  G  relative  to  the  bulk  modulus  K  was  found  from 
measurements  [7]  to  be  0.13.  In  Fig.  2  the  loss  tangent  in  G  is  assumed  zero. 

A  striking  feature  in  Fig.  2  is  the  appearance  of  a  spike  in  the  h/X  = 

0.5  curve  indicating  100/?  reflection.  This  may  be  explained  by  considering 

Fig.  3,  which  is  the  locus  of  points  for  which  the  relative  wave  speed  7  = 

c/c  is  such  that  zero  transmission  occurs,  as  a  function  of  the  dimensionless 
s 

frequency  k  d.  The  expression  for  this  was  given  in  Eq.  (11).  According  to 
s 

the  coincidence  relation  cq  =  c  sin  6,  the  points  for  7  when  the  angle  6  is 
varied  lie  on  a  vertical  line  starting  at  the  point  with  abscissa  k^d  and 
ordinate  c  /c  .  If  this  line  intersects  the  characteristic  curve  of  Fig.  3, 
total  reflection  will  occur.  In  the  example  depicted  in  Fig.  3  this  occurs 
for  7  =  18.9*  and  6  =  9.16*.  This  is  relatively  rare  for  near-pc  materials. 

A  small  amount  of  loss  in  the  shear  modulus  reduces  the  height  of  the  spike  to 
a  negligible  size. 
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A  global  viewing  of  the  nine  sets  of  curves  in  Fig.  2  shows  that  at  the 
intermediate  range  of  6,  say  from  30  to  60*,  there  is  little  difference  in  the 
appearance  of  the  curves  for  varying  density  and  sound  speed.  As  to  the  small 
angle  range,  it  appears  that  the  relevant  parameter  is  the  product  pc\  the 
groups  that  have  similar  appearance  are  a;  b,d;  c,e,g;  f,h;  and  i.  This 
corresponds  to  the  fact  that  for  fluids  the  behavior  near  0  =  0*  is  determined 
by  the  product  pc.  Near  6  -  90*  one  sees  that  there  are  three  groups 
distinguished  by  the  presence  or  absence  of  sero  reflection.  For  c0/c,j  = 
0.995  there  is  an  angle  where  the  reflection  is  sero;  for  CQ/C^  -  1-000  the 
sero  reflection  occurs  at  90*;  and  for  co/cj  =  1.005  there  is  no  angle  near 
90*  at  which  the  reflection  coefficient  is  sero. 

In  Figure  4  examples  are  given  where  both  the  density  and  the  sound  speed 
differ  by  + 5 %  or  - 5 %  compared  with  ideal  pc.  It  appears  that  too  small  a 
density  and  sound  speed  are  more  deleterious  than  too  large  a  density  and 
sound  speed,  judging  by  the  reflection  coefficient. 
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Fig.  4  -  Reflection  coefficient 
as  a  function  of 
incidence  angle  6 . 
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In  Figure  5  the  effect  of  a  finite  loss  tangent  in  G  is  shown,  namely 
=  0.01,  which  is  about  the  value  measured  for  the  shear  modulus  of  the 
rigid  pc  material.  Three  combinations  of  density  and  sound  speed  were  chosen 
One  sees  that  the  effect  is  quite  small. 
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